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Abstract 

m ; 

Final state interactions in the 5-wave hit system (1=0,2) are re-examined 
on the basis of the Omnes-Muskhelishvili equation and the coupled chan- 
Q^ ' nel formalism. The contributions to the pion scalar form factor from p 

and /2(1270) exchange in the t-channel and from the /o(980) s-channel 
Qh! resonance are separately evaluated and the role of the nontrivial polyno- 

mial in the Omnes function in a coupled channel situation is elucidated. 
Applications are made to K — * tttt and pp — > pp-inr. It is found that the 
contribution from the /o resonance to the form-factor is strongly reduced 
. ^ by a nearby zero. 

X' 
=3: 

1 Introduction 

Final state interaction (FSI) in the tot system plays an important role for many 
production reactions and meson decays. A case of long-standing interest is the 
AI = 1/2 rule in K — > trn decays. The experimental ratio of the decay amplitudes 
A j with isospin I — 0, 2 is [|T| 

- m) - 22 (1) 
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The calculated ratio is smaller [[| by at least a factor of 3 where this result 
includes perturbative QCD and soft-gluon corrections at the weak interaction 
vertex but no long-distance tttt FSI. In this paper we shall discuss the pionic FSI 
in the S- wave aiming at a concrete application to the AI = 1/2 rule for the 
K — > ittt decay and the pion production reaction pp — > Ttnpp. Our analysis shows 
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general features of FSI's which are relevant to other reactions involving pions or 
other hadrons. 

Several methods for the evaluation of FSI have been used in the literature. 
In one approach rescattering diagrams are evaluated directly. At low energies 
this has been done by applying chiral perturbation theory (CHPT) [||, |J]. The 
relevant application in our context is the calculation of the scalar form-factor of 
the pion in next to leading chiral order at low energies g |. To extend the 
calculations to s ~ 1 GeV 2 s-channel resonances and the coupling to the KK 
channel must be included. As a general tool the dispersive method based on the 
Omnes-Muskhelishvili (OM) equation f?|, || has turned out to be very efficient. 
It exploits analyticity and unitarity in order to connect the production or decay 
amplitude (or its form-factor) with the amplitude of elastic tttt scattering. To 
solve the OM equation we shall take the scattering phases either from phase shift 
analysis or from a theoretical model. We shall choose a model which satisfies 
the requirements of unitarity and analyticity, and hence the OM equation auto- 
matically. The model with parameters fitted to the experimental constraints is 
described in Sec. 3. 

For the K — > tttt decay it was realized a long time ago that the non-perturbative 
long-distance effects must be included. An enhancement of about a factor of 2 
in the 1 = amplitude was estimated to result from the broad cr(J PC = ++ ) 
meson ||. For K — > tttt the attraction in the 1 = channel must be combined 
with the repulsion in the 1 = 2 channel which favours the AI =1/2 rule. The 
analysis was done in CHPT to one loop in || |I0 . Rescattering in simple poten- 
tial model was evaluated in |TI|, [12| without regard to the energy dependence of 
the form-factor. An extensive study of the FSI effects in the S"-wave tttt system 
in production reactions and J/ip and if)' decays was conducted in ]T3|, [14], [15], CT . 



Unitarity and analyticity of the production amplitudes was taken into account 
in a self-consistent way. It was noticed, in particular, that a narrow resonance 
(/o in the present notation) in the tttt scattering phase 5j1q(s) corresponds to a 



shoulder in the tttt effective mass distribution in the reaction pp — > ppmr |rj, |T5]j . 
The occurrence of a shoulder rather than a peak results from an interplay of the 
resonant pole and a nearby zero. We shall discuss this feature in detail in Sec. 3. 
Resonance phenomena in the tttt S-wave were emphasized in ]17) where the fo 
resonance was discussed within a single-resonance model for the decay of a light 
higgs boson. The prediction of a drastic enhancement due to the fo resonance 
is in striking contrast with the findings for the pp — > ppmr reaction in |15 |. An 



analysis of the tttt final state interaction in the framework of the coupled channel 
OM equation was performed in U for the decay of a light higgs boson decay 
H —>■ 7T7T. In this evaluation the fo resonance also produced significant effects far 
below the KK threshold. 

The dynamics of the 1 = S-wave tttt interaction is characterized by several 
overlapping resonances [|l|, [16|, |I8| , narrow and broad. In the present paper 
we shall analyze the relative importance of the dynamical mechanisms in tttt 
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scattering for the calculation of the form-factors occurring in meson decays and 
in the pion pair production in pp scattering. In Sec. 2 we prepare the ground 
with an evaluation of the OM equation for a restricted energy range (the cut-off 
used excludes the /o resonance). With respect to the pion dynamics we shall 
mainly use the picture of |19[ which combines the p and f'2 exchanges in the 
t-channel with the fo resonance in the s-channel. The phases of the 7 = 0,2 
5-wave scattering are reproduced quite accurately in this model. To understand 
the role of the fo resonance for the calculation of the form-factor in the 1 = 
channel we shall introduce a coupled channel ansatz in Sec. 3. The final state 
interaction effects in the K — > nn decay are evaluated and the conclusions are 
presented in Sec. |. 

2 Form- factors from the Omnes — Muskhelishvili 
equation 

The OM equation |7|, |§ connects the form-factor F(s) with the elastic final state 
scattering phase 5(s). For a single channel problem the OM equation is 

where a once-subtracted form has been used. The general solution of (Q) has the 
form 



(3) 



as long as 



\F(s)\ 

6(s) — > const, - — — ^ for s — > 00. (4) 

The polynomial P(s) is real for s real. In special cases, like potential scattering 
without bound states, P(s) is a constant, but in general additional information 
is required to determine it. 

For K — > 7T7T in the simplest evaluations single channel tttt scattering data are 
used below the KK threshold (the coupling to the 4n channel is known to be 
small). Fig. [l|a shows the mi J = I = scattering phase Sq(s) from the phase 
shift analysis ||15|| . In Fig. |l|b we show the same phase from the meson exchange 
model mentioned earlier and developed in |19|, |2T| . We briefly recapitulate 
the ingredients for the benefit of the later discussion. The phases in Fig. [TJb 
correspond to unitarized p and /2 exchange with the fo resonance added. The 
individual contributions are shown in the figure as explained in the capture. The 
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Born term for the p-exchange is 



ns,ty B = A 



2G 



u 
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s-t 



t 



u 



(5) 
(6) 



where m p is the mass of the p meson, G = g> 2 /327T, and g pwn is the pirn coupling 
constant. The 1 = 2 amplitude will be needed later. The S'-wave projection is 

- 4m!" 



Tba°-s( s ) 
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(7) 



X-matrix unitarization is introduced by 



Kl(s) 



where 



Tba-s( s ) 



(8) 



(9) 



and p(s) = (1 — 4m 2 /s) 1 / 2 . The coupling constant g p7r7r is determined from the 
p meson decay width in the crossed 1 = 1 channel after i^-matrix unitarization 
nj. The corresponding value is g p7T7r = 6.04 which is close to the result obtained 



from the KSFR relation 25 



The corresponding expression for f 2 exchange in Born approximation |L9 

-11 2 
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2G 



fit 
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(2s + m\ - Ami) 2 ~ ( m % - 4m 2 ) 2 /3 
s — Ami 



ln(l + 



4m 2 



(10) 
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where G h ~ 0.19 GeV~ 2 . 

The /o resonance is included using the Dalitz-Tuan representation [ffiH , i- e - 
the S-matrix is considered to be the product of the S-matrices corresponding to 
the individual mechanisms. The corresponding Breit-Wigner parametrization is 



taken from 19 



S(s) 



s — M 2 — ig 1 p 1 (s) + ig 2 p 2 {s) 
s - M 2 + ipipi(s) + ig 2 p 2 (s) 



(12) 



where pi(s) = yl — 4m 2 /s, p 2 (s) = yl — Am 2 K /s, and the resonance parameters 
are M r = 0.9535 GeV, g x = 0.1108 GeV 2 , g 2 = 0.4229 GeV 2 , respectively. The 
scattering phase in the meson exchange model is a good description of the data 
for s < 1.4 GeV 2 , see Fig. [TJb. 
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Figure 1: The S-wave scattering phase 5q vs. ^Js: (a) the Kl fit from [ |15|] , 
(b) the meson exchange model described in the text (solid line: the total phase; 
dashed line: P + /2 ^-exchange; dotted line: fo resonance). The experimental data 
are from [E2L 



In order to show the sensitivity of the form-factors to the variations in the 
phase we shall evaluate the OM equation using the two sets shown in Fig. |l]a and 
|l|b. For small s the integral of Eq. ([$]) is dominated by low energies, see 0. As a 
first step we evaluate in this section the OM equation with a cut-off A. Choosing 
A = 0.975 GeV we exclude the fo and the KK threshold region as in some early 
applications. We therefore write 

F{s) = F A (s)P(s) (13) 

where 

Fa(s) = exp U f -^Pr^ds) . (14) 

\7T J Ami S'{S' — S) J 

The polynomial P(s) represents the contribution from high energies and any 
other dynamics not included so far. We observe that any additive contributions 
in the phase lead to a multiplicative factor in F, see eq. ([3|). We shall use 

P(s) = l + bs (15) 

where the parameter b is related to the scalar radius of the pion by 

F(s) = l + ±(r 2 s )s. (16) 
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Figure 2: The pion scalar form-factor F I=0 (s) vs. y/s evaluated using the OM 
equation with the cut-off A = 0.975 Gev and the mr scattering phases as shown 
in Fig. |l]a and [l|b, correspondingly. 



When plotting Fig. ^| we have adjusted the polynomial (]i~5[) in order to have 
(r 2 ) = 0.6 fm 2 in both cases. This leads to b = 0.32 GeV -2 for Fig. |a 
and b = 0.83 GeV -2 for Fig. ^jp. Figure shows that the form-factor is rather 
sensitive to the irn scattering phase shift, one of the important parameters being 
the scattering length o!q. For the meson exchange model a[j = 0.24m" 1 and the 
absolute value of the form-factor continues to rise between the tttt threshold and 
\/s ~ 0.5 GeV. On the other hand the form-factor corresponding to the phase in 
Fig. |2|a displays a more prominent cusp due to a larger value of a[j = 0.51m" 1 
and decreases above the -rm threshold. For y/s < 0.5 GeV the result in Fig. ||b is 
close to the solution of the coupled- channel OM equation in ||. 

3 Protective Zero and the fo Resonance 

For energies around s = 1 GeV 2 the truncation in the calculation of the form- 
factor must be abandoned and the role of the /o and the KK threshold discussed. 
The resonant part of the phase will be defined by 



Arctan 



gk(s) 
(M 2 - s) 



(17) 
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corresponding to the resonance amplitude 



gk(s) 

■L res (S ) 



s - M 2 + igk(s) 



where k(s) = Js — 4m 2 /2, leads to 




(V ' J ' S,) ds') = (19) 



7T Js S'(s' — s) 

M 2 + gm^ 



M 2 - s - igk(s) 

Inserting this phase naively into Eq. ( fTi]) has the undesirable feature that 



(20) 



\F res (s)\ S -°° (21) 

rather than unity which would be expected at high energies where the resonance 
contribution should vanish. The wrong asymptotic form is actually imposed on 
the whole solution by means of phase additivity. In Fig. |3| the dotted line corre- 
sponds to the naive evaluation of F(s), with P(s) being set to unity. Apart from 
the wrong asymptotics it is also seen that the fo resonance dominates the form- 
factor far outside the resonance region Mf ± Tf . Recall that the experimental 
width is Tf « 60 MeV. It is clear that this defect should be compensated by a 
non-trivial polynomial P(s) in the solution of the OM equation. 



3.1 The OM equation for a resonance in the Weisskopf- 
Wigner model 

We shall study the modification required for a sensible inclusion of a direct chan- 
nel resonance into the OM equation by means of a very simple coupled channel 
model. The following nonrelativistic ansatz, which is a variant of the Weisskopf- 
Wigner (WW) model, already has all the necessary ingredients. Only one scat- 
tering channel is introduced which has no diagonal potential. It will be denoted 
by its relative momentum \k). The only interaction in the model results from 
the coupling to a bound state \b) (effectively representing a second channel). We 
assume a channel coupling of the form 

(k\V\b) = rf(k) = -^— 2 (22) 

where 7 is the coupling constant (dimension [7] = [A;] 3 / 2 ) and /x characterizes the 
range of interaction. The T-matrix satisfies the Lippmann-Schwinger equation 

T{E) = V ^ V{\ + G (E)T(E)) (23) 
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Figure 3: The pion scalar form-factor F I=0 (s) vs. \/s. The solid line corresponds 
to p and f'2 exchange and an /o resonance including the polynomial (protective 
zero). The f resonance alone leads to the dotted line (OM equation without a 
polynomial) and to the dashed line with polynomial. 



where Gq(E) is the free Green function and E = k 2 /2m (m is the reduced mass). 
The solution for the scattering amplitude has the form 

f(k) = -2m(k\T{E)\k) = (24) 
-2m 7 2 £ 2 (A;) 



^-E b -l 2 D(k) 



(25) 



with 



771 

D{k) = (b\VG (E)V\b) = -jfi^Tjy ■ (26) 

In our model the form-factor F(k) describing the final state interaction is 
equal to the scattering wave function at zero distance according to standard 
results from scattering theory Wt 



F(k) = (r = 0\k {+) ) = (r = 0\k) + (r = 0\G (E)T(E)\k) (27) 

-f 2 Z(k)£(k) . . 

= 1+ p USU (28) 

±--E r -^Dik) 
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with 



Z{k) = (29) 

In Fig. |] we show the scattering phase and the form-factor for the WW model 
for fj, = 5m, 7 = 10/z 3//2 , E r = 4m. The pole produces a resonance peak in the 
energy dependence of the form-factor which is damped by a nearby zero restoring 
the right limit F — > 1 for E —>■ oo. The reduction imposed by this zero is 
enormous. 




Figure 4: The momentum dependence of the scattering phase (a) and the form- 
factor (b) for a resonance in the WW model. The solid line is the exact solution 
(|28|), the dashed line is the solution of the OM equation without polynomial 
factor. The dots show the approximate solution of the OM equation with the 
factor (E — E z ). 



Formula ( p8|) can be rewritten explicitly showing the interplay of the pole and 
the zero: 

n*> = HI (30) 

^y 2 m I 

A{E) = E — E r — — (31) 

fx (k 2 + fx 2 ) 

Rm F F I 2 ™ (k 2 -ii 2 - 2ink) 

B{E) ~ E - Er ~ll — WT^y ■ (32) 
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In the limit of weak coupling the resonance in the scattering channel is directly 
connected to the bound state in the continuum which has an energy shift AE r 
and a width T r : 

7 2 m(k 2 - fi 2 ) 

r = ^ mk (34) 

(k 2 + y 2 ) 2 [6V 

where E r = k 2 /2m. The form-factor in the vicinity of the resonance has the form 



F{k) = E-(E r + AEr-iT r /2) (35) 
and the zero E z is located near the resonance energy E T + AE r at 

E z = E r + AE r + j-T r (36) 

If \E Z — E r — AE r \ > T r , the resonance produces a pronounced peak followed by 
a dip in the energy dependence of the form-factor. In case \E Z — E r — AE r \ < T r 
the energy dependence coming from the pole is damped completely by the zero in 
the nominator, and only a dip is visible in the form-factor. Notice that the zero 
is of dynamical nature and disappears for vanishing channel coupling: F — > 1 as 
7^0. 

Since A(E) is a real symmetric function of momentum k, it does not contribute 
to the elastic scattering amplitude. The solution of the OM equation without a 



polynomial factor reflects only the resonance pole in formula fl35|) as shown in 
Fig. £|, dashed line. By including the factor (E — E z ) one gets 



E z \ir Jo E'(E' - E) 



which is very close to the exact solutionQof the WW model. 

These results characterize a coupled- channel resonance. The scattering phase 
beyond the resonance does not decrease as it occurs for a direct channel poten- 
tial resonance^] where no extra polynomial factor appears in the solution of the 
OM equation. For a potential resonance the decrease of the phase for s —>■ oo 
guarantees that the asymptotic limit of the form-factor is one. 



1 A careful analysis of the OM equation for the model considered shows that there is an extra 
factor (fc 2 + v 2 )/(k 2 + fi 2 ) resulting from the singularities in the upper halfplane of complex 
momentum k: a pole at k — i[i and a nearby zero at k — iv. For our example this factor is 
close to 1 in the region of the resonance. 

2 In the literature the first category is often called normal resonance and the second one 
molecular or bootstrap resonance, see e.g. J15| and references therein. 
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It must be emphasized that in the WW model considered, the resonance-dip 
structure occurs only in processes where the particles in the scattering channel 
are produced at small distance due to some extraneous interaction which can 
be treated perturbatively, so that the momentum dependence of the production 
amplitude is entirely determined by the form factor F(k) given by Eq. ( |28| ) (this 
is relevant for the K —>■ tttc decay). This situation must be distinguished from 
a situation where the original bound state \b) is produced as a resonance with 
amplitude C and then decays into the scattering channel. The corresponding 
amplitude with rescattering included is 

T b (k) = C— 2 ^ (38) 

2m 

which has a purely resonant behavior, there is no nearby zero. Studying the 
energy dependence of the data in the vicinity of the resonance one can determine 
whether this situation is realized for the process in question. 



3.2 Application to the /o resonance and constraint from 

PP — > pp7T7T 

To evaluate the role of the fo resonance for K — > tttt decay we use the S-matrix 
in Breit-Wigner form fitted to data |19[], see Eq. ([12]). As we demonstrated in 
Sec. PTT] , the polynomial in the solution of the OM equation is expected to have 
a zero at s = s z close to the resonance: 

P(s) = 1 - - (39) 

Note that s z — > as g± — > 0. In order to fix the position of the zero we 
use information from a related process and study the effective mass distribution 
(M = -y/i) of pion pairs produced in the reaction pp pp-rni p8 |, which can be 
expressed by fl4[ 



*L ~ {M2 -^ in \F<M>)? (40) 

Including the polynomial ([39]) into the calculation of the form-factor F(s) (fo 
plus p and f% exchange) we obtain s z = 1.0 GeV 2 for the position of the zero, 
see Fig. |5[ The fit shown for the mass distribution da/dM also contains a factor 
(1 + 0. 25s) in the polynomial and an overall normalization constant. The position 
of the zero however, is determined very precisely from nearby data alone. The 
corresponding scalar form-factor will be discussed in sec. [|. 

Note that the factor containing the zero can be incorporated into a formal 
solution of the OM equation, if a physically equivalent discontinuous scattering 
phase is introduced 

5( s ) = 5(s) - tt0(s - s z ) . (41) 
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Figure 5: The effective mass distribution of pion pairs in pp — > ppim vs. M 



^Js. The data are from |28 



It is interesting to note that the description of the pion pair distribution da / dM 
in [14] , llEf was seemingly achieved using the trivial polynomial P(s) = 1. However 
the elastic phase was calculated from the expression $ = Arctan(Im Tn/Re Tn). 
In the presence of inelasticities the phase of T n is bounded to the interval [0, tt) by 
the requirement of continuity. When Re T n changes from negative to positive due 
to the sharp resonance rise of 5(s) the phase $ drops sharply by nearly n. With 
this choice the Omnes function develops a zero close to the point where 5 = tt. 
While this provides a good description of the data, the introduction of the zero 
in this way appears to be accidental. For instance, in the model considered in 
Sec. |3.1| there is no connection between the position of the zero and the condition 
5 = 7i. Also, if the scattering phase 6 reached tt before the KK threshold, the 
zero factor would not be obtained from using the phase prescription for $ quoted 
above. 



4 Summary of results for the K —> nn decay and 
conclusion 

For the AI — |, | K — > 2tt decays, without long-range final state interactions, 
the bare weak vertex can be parametrized as 0, T b I=02 {s) = Ci=Q^{.s — m^) where 
s = m? K for an on-shell kaon. The quantity T b contributes to the polynomial in 
the Omnes representation of the K — > ixix decay. We have prepared the ground for 
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the 1 = S-w&ve final state interaction in the preceeding sections. The solid line 
in Fig. |^ shows the net result for the model combining the p and /2 exchange with 
the fo resonance. The resonance in the form-factor is protected by the zero at 
s z = 1 GeV 2 as determined from the pion pair production data. At the kaon mass 
the 1 = enhancement factor is F{m 2 K ) = 1.62, a result which is similar to the 
values obtained in the literature quoted above. From p exchange alone we obtain 
F(m 2 K ) = 1.38, p and fi give F{m 2 K ) = 1.57 while the enhancement from the fo 
resonance alone is F(m 2 K ) = 1.03. For the complete form-factor the reduction of 
the fo contribution induced by the protective zero is of course crucial. The effects 
of the zero and the resonance largely cancel and only a very small contribution to 
the form-factor far away from the pole (zero) survives. For example, at s = the 
pion scalar radius is: (r 2 ) = 0.52 fm 2 when only considering p and f2 exchanges. 
When including the resonance protected by the zero we have (r 2 ) = 0.58 fm 2 . 
We see that the inclusion of the resonance does improve the result on the scalar 
radius but avoids too large an effect. Our full result is very close to the value 
obtained in M where (r 2 ) is determined from chiral perturbation theory. Without 
the zero we would have obtained a rather large value (r 2 ) = 0.81 fm 2 . 

In order to complete the evaluation of the the overall AI = 1/2 enhancement 
factor the contribution of the 1 = 2 channel must be evaluated as well. Due to 
the signature of the crossing matrix the contribution from p exchange is repulsive 
in the 1 = 2 channel, see (^[). On the other hand fa exchange does not change 
sign relative to the 1 = channel leading to destructive interference between 
p and /2 for the isotensor. The solid line in Fig. ^]a shows the unitarized sum 
of p and /2 exchange. Also shown is p exchange modified by a vertex form 
factor with monopole range A p = 1.5 GeV (dashed line) which is a good effective 
parametrization of the data. The phases at higher energies are not known, but 
fortunately the form-factor at y/s = rriK is not sensitive to this region. The 
corresponding form-factor F I=2 (s) is shown in Fig. |6|b. At the kaon mass we 
obtain a reduction factor F I=2 (m 2 K ) = 0.9 leading to a combined AI = 1/2 
enhancement of F I=0 {m 2 K ) / F I=2 {m 2 K ) = 1.81 which is satisfactory, but slightly 
less than the value required by the data. 

We conclude that the p and /2 exchange interactions remain the dominant 
mechanisms for the FSI enhancement factor in the AI =1/2 rule in K — > im. 
Since p exchange generates a broad pole in the 1 = S-wave amplitude [jl9| 
one can associate this enhancement with a a meson. The fo resonance plays a 
minor role. This is due to the occurrence of a protective zero at s = 1 GeV 2 
modifying the polynomial in the OM equation. The nature and position of this 
zero has been verified by analyzing pion pair production in pp ppmr where the 
fo resonance only leads to a small shoulder in the mass distribution. A simple 
coupled channel model describes this situation very adequately. We expect that 
this damping mechanism will be applicable to many other decay and production 
reactions in the vicinity of a coupled channel resonance. 
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Figure 6: The tttt I = 2 S-wave scattering phase 5q 2 vs. y/s (a) and the form- 
factor F I=2 (s) (b). Solid line: p + f'2 exchange, dashed line: p-exchange with 



vertex form-factor. The experimental data are from |p9 | 
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